I Applications 
in Plant Sciences 



Applications in Plant Sciences 2014 2(9): 1400047 

Primer Note 



mlcrosatellite primers in agave utahensis 
(Asparagaceae), a keystone species in the Mojave 
Desert and Colorado Plateau 1 

Charlee Byers 2 , Peter J. Maughan 2 , Jared Clouse 2 , and J. Ryan Stewart 2 ' 3 

department of Plant and Wildlife Sciences, Brigham Young University, 4105 Life Sciences Building, Provo, Utah 84602 USA 

• Premise of the study: Utah agave (Agave utahensis) and its putative subspecies, A. utahensis subsp. kaibabensis and A. utahen- 
sis subsp. utahensis, are keystone species of the Mojave Desert and Colorado Plateau in the southwestern United States. Here 
we developed microsatellite markers to study population structure and genetic diversity of the two subspecies of A. utahensis. 

• Methods and Results: We analyzed 22,386 454-pyrosequencing large contigs (>400 bp), derived from a genome reduction 
experiment consisting of A. utahensis accessions, for putative micro satellites. The use of unique multiplex barcodes for 
each of the Agave accessions allowed for the identification of putatively polymorphic microsatellites based solely on 
sequence alignment analysis. We report the characteristics of 1 1 polymorphic microsatellite loci based on a panel of 104 
individuals from the two subspecies. The number of alleles per locus varied from three to eight, with an average of 5.5 
alleles per locus. Observed and expected heterozygosity values ranged from 0.038 to 0.777 and 0.038 to 0.707, 
respectively. 

• Conclusions: The microsatellites identified here will be invaluable for future studies of population structure, polyploidy, and 
genetic diversity across the species. 

Key words: Agave utahensis subsp. kaibabensis; Agave utahensis subsp. utahensis; Asparagaceae; microsatellite; simple 
sequence repeat (SSR). 



Agave utahensis Engelm. (Utah agave; Asparagaceae) and its 
putative subspecies, A. utahensis subsp. kaibabensis (McKelvey) 
Gentry and A. utahensis subsp. utahensis, play critical roles 
as keystone taxa in their native habitats. Utah agave contrib- 
utes significantly to soil formation and the creation of oases of 
organic matter in the Mojave Desert and Colorado Plateau in 
the southwestern United States (Gentry, 1982). Their large in- 
florescences produce copious amounts of nectar for a wide array 
of pollinators, including bats, birds, bees, and hawkmoths 
(Slauson, 2001). The rise in frequent and intense fires across 
the distribution of Utah agave is threatening native popula- 
tions and consequently the genetic diversity of the species. As 
with other agaves, Utah agave reproduces clonally via vegeta- 
tive rhizomes (Gentry, 1982; Arizaga and Ezcurra, 2002). As 
such, the risk exists that its asexual reproductive strategy 
could hasten its ecological decline if fire frequencies continue 
to increase. 

Moreover, taxonomic questions exist concerning the sub- 
species ranking within this highly variable species complex. 
The subspecies ranking of the two Utah agave (subsp. kaiba- 
bensis and subsp. utahensis) is based solely on morphologi- 
cal traits, with subsp. kaibabensis having wider and greener 
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leaves (3.0-5.5 cm) and larger rosettes (40-100 cm) than 
subsp. utahensis (1.0-3.0 cm, 15-40 cm, respectively) (Gentry, 
1982; Reveal and Hodgson, 2003; Janeba, 2008). DNA marker 
information will help characterize the genetic diversity of the 
Utah agave species complex and confirm whether subspecies 
ranking is justified. 

Despite its importance as a keystone species, no molecular 
tools have been developed to manage Utah agave or to inves- 
tigate population structure, subspecies ranking, hybridization, 
or polyploidy within the species. We constructed primers to 
amplify microsatellite markers from both subspecies along a 
west-east transect across southwestern Utah and north-central 
Arizona. 

METHODS AND RESULTS 

The genomic reduction protocol used here is described in detail by 
Maughan et al. (2009) for single-nucleotide polymorphism discovery. Here 
we use the same methodology, but analyze the results for putatively poly- 
morphic microsatellite sequences. The method is based on the conservation 
of restriction sites across individuals followed by genomic selection and 
en masse fragment sequencing via next-generation sequencing technology. 
Specific DNA barcodes are incorporated into the DNA fragment of each 
member of the screening panel, which are subsequently used to deconvolute 
the sequence read pool and to search for polymorphic microsatellite loci 
across the assembled DNA fragments. Leaf tissue was lyophilized and DNA 
extracted using a modified mini-salts protocol as reported by Todd and Vodkin 
(1996). In short, 450 ng of DNA from two A. utahensis subsp. kaibabensis 
and two A. utahensis subsp. utahensis agave accessions were separately 
double-digested for 1 h at 37°C with 3 units of the restriction enzymes 
EcoRl and Bfal (New England Biolabs, Beverly, Massachusetts, USA) in 
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lx NEB4 restriction buffer. The resultant DNA fragments were immediately 
ligated with 1.5 uM 5'-TEG biotinylated/3'-phosphorylated EcoRl adapters 
and 15 jiM S'-phosphorylated Bfal adapters (see Maughan et al., 2010 for 
adapter sequences) using 3 units of T 4 ligase (New England Biolabs) at 
16°C for 3 h. Non-biotin-labeled DNA fragments (Bfal to Bfal restriction 
fragments) were removed from the reaction using a biotin-streptavidin para- 
magnetic bead separation using M-280 streptavidin beads (Invitrogen, 
Carlsbad, California, USA) according to the manufacturer's specifications. 
The remaining EcoRl-Bfal and EcoRl-EcoRl fragments, with the biotin la- 
bel still attached, were resuspended in 100 jiL of 10 : 1 TE buffer (10 mM 
Tris, 1 mM EDTA, pH 8.0). 

Four primer pairs, designed to be complementary to the EcoRl and Bfal 
adapter sequences and to carry unique 5' 10-base barcode sequences 
(Maughan et al., 2010), were synthesized by Integrated DNA Technologies 
(Iowa City, Iowa, USA). A single barcode primer pair was used to amplify 
1 JiL of the streptavidin-cleaned DNA fragments for each DNA sample. 
Amplification of each sample was performed in 50-uX PCR reactions using 
lx Advantage HF 2 PCR Master Mix (ClonTech, Mountain View, Califor- 
nia, USA) and 0.2 uM of the MIDX-£coRI and MIDX-Bfal primer pairs, 
with the following thermocycling profile: 95°C for 1 min followed by 22 
cycles of 95°C for 15 s, 65°C for 30 s, and 68°C for 2 min. Equimolar 
amounts of DNA for each of the PCR reactions were pooled and sequenced 
using standard protocols for 454-pyrosequencing as a service at the Brigham 
Young University DNA Sequencing Center (DNASC, Provo, Utah, USA) 
using a Roche-454 GS FLX instrument and Titanium reagents (454 Life 
Sciences, a Roche Company, Branford, Connecticut, USA) without DNA 
fragmentation. DNA reads were bioinformatically trimmed and separated 
into barcode pools representing the four Utah agave accessions and then de 
novo assembled together using the default parameters of Roche Newbler as- 
sembler (version 2.3; 454 Life Sciences, a Roche Company). All 786,104 se- 
quence reads were assembled into 22,386 large reference contigs (>400 bp). 
Using the computer program MISA (Thiel et al., 2003), the reference contigs 
were scanned for perfect di-, tri-, and tetranucleotide motifs with minimum 
repeat units of eight, six, and five, respectively . 

A total of 1013 putative microsatellite markers were identified, including 
762 di-, 205 tri-, and 46 tetranucleotide motifs. The most common motifs 
were AG/CT, AAG/CTT, and AAAC/GTTT. The MISA output files and the 
reference contig.ace file were then processed using custom Perl scripts 
(Stajich et al., 2002) to preselect microsatellite loci that varied in repeat 
length as determined from consensus reads for each of the agave subspecies (uta- 
hensis vs. kaibabensis). A total of 298 putatively polymorphic microsatellites 



were identified, of which 15 trinucleotides representing the range of motifs 
(AAT, ACA, CAC, CCA, GAG, GCT, GTT, TAT, TTC) and repeat lengths 
(5-10 repeats) were selected for primer design using Primer3 version 2.0 
software (Rozen and Skaletsky, 2000) and validation on a diversity panel 
of 104 Utah agave accessions collected from eight populations in southern 
Utah and northern Arizona (Appendix 1). Amplification of each sample 
was performed using multiplex PCR reactions, which included three prim- 
ers, each with different fluorescent tags (VIC, 6-FAM, NED; Table 1). Mul- 
tiplex PCR reactions were performed using a Type-it Microsatellite PCR Kit 
(QIAGEN, Germantown, Maryland, USA) in 10-jiL reactions with 30 ng 
of genomic DNA and 0.2 jiM of all primers according to the manufacturer 
guidelines with the following amplification parameters: 95°C for 5 min, fol- 
lowed by 28 cycles of 95°C for 30 s, 57°C for 90 s, 72°C for 30 s, followed 
by a final extension cycle of 60°C for 30 min. The amplified products were 
diluted 1 : 60, and 1 uL of the dried products were analyzed using frag- 
ment analysis on an ABI 3730x1 using GeneScan 500 ROX (Applied Bio- 
systems, Carlsbad, California, USA) as the size standard at the Brigham 
Young University DNASC. The data generated from the fragment analysis 
were scored using the microsatellite plugin in Geneious R6 (Drummond 
et al., 201 1) and analyzed using default parameters of the computer program 
Arlequin version 3.5 (Excoffier and Lischer, 2010). 

Four of the 15 preselected markers failed to amplify or produced poor, weak, 
or ambiguous amplification products. Eleven markers were polymorphic, pro- 
ducing a total of 61 alleles across both subspecies. The number of alleles per 
locus varied from a low of three alleles (BYU5866) to a high of eight alleles 
(B YU5 164 and BYU8677), with an average of 5.5 alleles per marker (Table 2). 
Interestingly, two of the microsatellite markers (BYU4463 and BYU4988) ex- 
hibited individuals (-25%) that amplified between two and four alleles, pre- 
sumably due to the amplification of two paralogous DNA sequences. The mean 
observed and expected heterozygosity values for A. utahensis subsp. utahensis 
and A. utahensis subsp. kaibabensis were 0.357 and 0.491 and 0.357 and 0.408, 
respectively (Table 2). Division of the diversity panel by putative subspecies 
ranking exhibited a between-subspecies genie differentiation (F ST ) of 0.24. 
When all individuals were analyzed for Hardy- Weinberg equilibrium (HWE) 
on a subspecies basis, irrespective of subpopulation, several markers violated 
HWE (Table 2). This was likely due to population substructure (i.e., Wahlund 
effect). Indeed, when the data were analyzed for each of the eight subpopula- 
tions (Appendix 1), all microsatellites were in HWE for five of the populations, 
with only BYU4012 or BYU8677 violating HWE in the remaining three popu- 
lations (note that the sample sizes for these populations are too small for defini- 
tive conclusions). 



Table 1. Characteristics of 1 1 polymorphic microsatellite primers developed for Agave utahensis. 







Allele size Multiplex 


BLASTX similarity 


GenBank 


Locus 


Primer sequences (5'-3') 


Repeat motif range (bp) T a (°C) set 


(£- value) a 


accession no. 



BYU3268 


F 
R 


VIC-CCTCCAATCAAATTCCTCA 
TGAAAGACCCCTCATCATCC 


(GCT) 5 


138- 


-147 


57 


2 


IAA-leucine resistant-like gene 
XP_007044503.1 (2e-29) 


KJ855288 


BYU3674 


F 
R 


VIC-CATCAATTAATAACATGAAGAATGG 
TGGTTCATGTAGCCTACCCC 


(ACA) 6 


147- 


-156 


57 


3 


Elongation factor 1 -alpha 
XP_003581610.1 (le-24) 


KJ855289 


BYU4012 


F 
R 


VIC-TTTTTGGGAGTAAGCTGGGA 
CCCGGGGAGTAAAGAGAGAG 


(CAC) 5 


131- 


-140 


57 


1 


N/A 


KJ855290 


BYU4463 


F 
R 


VIC-CGTTGAGCTTGTTGAAGTCG 
AACCACACCAAAGTCAAGCAG 


(GAG) 6 


178- 


-199 


57 


4 


Allene oxide synthase 1 
NP_001 105244.1 (le-68) 


KJ855291 


BYU4988 


F 
R 


NED-GACGCTGGTATATTTCTCCTG 
TCTCGGATCTTCAATCTGGG 


(TAT) 5 


172- 


-190 


57 


1 


N/A 


KJ855292 


BYU5164 


F 

R 


6FAM-ATCGCTCCTCCAGAGTTTCA 
GTGTGTGCTTTGTTGGTTGG 


(CCA) 7 


155- 


-176 


57 


1 


N/A 


KJ855293 


BYU5866 


F 

R 


6 FAM- ACACAGGCCCAACACAAAAT 
GGGACAATCCTTCCATTAAACA 


(AAT) 5 


164- 


-170 


57 


3 


N/A 


KJ855294 


BYU7269 


F 

R 


6FAM-GATCCATAGCCGACTCCAAA 
GAGGATGACAAACAATCCGC 


(GTT) 6 


155- 


-167 


57 


2 


N/A 


KJ855295 


BYU8490 


F 

R 


6FAM-CCAAATACTTGGGACAAATGG 
CAGAAATAGAAGGCCTTAATAAGTGG 


(GTT) 6 


173- 


-191 


57 


4 


N/A 


KJ855296 


BYU8677 


F 

R 


NE D-AAAC C AAAAAC AGGC C AC T G 
TGAGTCGAGGGAATTAGGGA 


(AAT) 10 


171- 


-198 


57 


2 


N/A 


KJ855297 


BYU9897 


F 

R 


NED-TAGTTAATTGCGCCCCCTTT 
CCGGAGTTGGGAGGAGTC 


(TTC) 6 


177- 


-189 


57 


3 


N/A 


KJ855298 



Note: N/A = no significant associated hits to the database; T a = annealing temperature. 
a BLASTX search against RefSeq database (£- value < le-15). 
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Table 2. Genetic properties, including number of alleles detected at 
each micro satellite marker and observed and expected heterozygosity 
values, for each subspecies. a 



A. utahensis subsp. A. utahensis subsp. 

utahensis (n = 78) c kaibabensis (n = 26) c 



Locus b 


A 


H 0 


H e 


A 


H 0 


He 


A 


BYU3268 


4 


0.3333 


0.3409 


4 


0.0384 


0.0384 


2 


BYU3674 


4 


0.3866 


0.3863 


4 


0.3043 


0.2695 


3 


BYU4012 


4 


0.1621 


0.3546 


4 


0.1600 


0.5355 


3 


BYU4463 


7 


NT 


NT 


7 


NT 


NT 


4 


BYU4988 


7 


NT 


NT 


6 


NT 


NT 


5 


BYU5164 


8 


0.5223 


0.6615 


8 


0.7777 


0.7031 


4 


BYU5866 


3 


0.2343 


0.4769 


3 


0.3913 


0.4492 


2 


BYU7269 


5 


0.5384 


0.6832 


5 


0.4230 


0.5316 


3 


BYU8490 


7 


0.3000 


0.4642 


6 


0.0416 


0.0416 


2 


BYU8677 


8 


0.4000 


0.7067 


7 


0.3846 


0.5399 


4 


BYU9897 


4 


0.3333 


0.3409 


1 


0.6956 


0.5642 


3 


Average 


5.5 


0.3567 


0.4905 


5.0 


0.3573 


0.4081 


3.2 



Note: A = number of alleles; F ST = differentiation between subspecies; 
H e = expected heterozygosity; H Q = observed heterozygosity; n = number 
of individuals sampled; NT = not tested. 

a Bold values indicate a significant deviation (P < 0.001) from Hardy- 
Weinberg equilibrium when analyzed by subspecies. 

b Microsatellite markers BYU4463 and BYU4988 were not tested for H Q 
and H e because more than two alleles were identified in several of the 
individuals analyzed. 

C F ST = 0.2407 (P< 0.0001). 

CONCLUSIONS 

We used a novel genomic reduction approach tied to 454- 
pyrosequencing and DNA sample barcoding to identify polymor- 
phic marker loci in A. utahensis. Eleven microsatellite markers 
were validated on 104 accessions of A. utahensis subsp. utahen- 
sis and subsp. kaibabensis. These markers should prove invalu- 
able for performing larger investigations of population structure, 
genetic diversity, and phylogenetic relationships in agaves. 
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Appendix 1. Geographic origin of the 104 Agave utahensis accessions studied. Plants sampled within a location were separated by a minimum of 10 m. 
Accessions were deposited in the Stanley L. Welsh Herbarium (BRY), Brigham Young University. 



Population location 


Subspecies 


N 


Accession no. 


Latitude 


Longitude 


Silver Reef 


utahensis 


14 


624001 


37°16 , 09.7 ,, N 


113 0 22'41.3"W 


Castle Cliff 


utahensis 


15 


624005 


37°03 , 55.0 ,, N 


113°53T6.0"W 


Woodbury 


utahensis 


17 


624009 


37°01'44.1"N 


113°52'11.7"W 


Lime Kiln Canyon 


utahensis 


15 


624011 


36°38 , 41.4 ,, N 


114°01'05.5"W 


Red Rock Spring 


utahensis 


6 


624016 


36°33 , 06.6 ,, N 


114°01'30.2"W 


Ledges Parkway 


utahensis 


11 


624023 


37°12 , 24.7 ,, N 


113°35'59.8"W 


Navajo Nation GC 


kaibabensis 


16 


624017 


35°56'14.1"N 


111°41'28.4"W 


Navajo Nation-Hoggan 


kaibabensis 


10 


624020 


35°55'42.2"N 


111 0 34'44.7"W 



Note: N= number of individuals sampled at each location. 
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